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^ ■ Abstract 

We show that expanding or contracting Kasner universes are unstable due to the amplifica- 
tion of gravitational waves (GW). As an application of this general relativity effect, we consider 
a pre-inflationary anisotropic geometry characterized by a Kasner-like expansion, which is driven 



. dynamically towards inflation by a scalar field. We investigate the evolution of linear metric fluc- 

Q \ tuations around this background, and calculate the amplification of the long- wavelength GW of 

H ■ 

' a certain polarization during the anisotropic expansion (this effect is absent for another GW po- 

: 

larization, and for scalar fluctuations). These GW are superimposed to the usual tensor modes of 
>^ , quantum origin from inflation, and are potentially observable if the total number of inflationary 

\ e-folds exceeds the minimum required to homogenize the observable universe only by a small mar- 

m ■ 

gin. Their contribution to the temperature anisotropy angular power spectrum decreases with the 
multipole d. as where p depends on the slope of the initial GW power-spectrum. Constraints 



o 

OO . 

, on the long-wavelength GW can be translated into limits on the total duration of inflation and 



the initial GW amplitude. The instability of classical GW (and zero-vacuum fluctuations of gravi- 
tons) during Kasner-like expansion (or contraction) may have other interesting applications. In 
particular, if GW become non-linear, they can significantly alter the geometry before the onset of 
inflation. 
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I. INTRODUCTION: ANISOTROPIC PRE-INFLATION 



The inflationary stage of tlie very early universe explains the dynamical origin of the 
observed isotropic and homogeneous FRW geometry. The patch of the FRW geometry 
covers the cosmological horizon and beyond if inflation lasted 

e-folds or longer. Here V is the potential energy of the inflation, and A is a correction 
from the (p)reheating stage after inflation, which is not essential for our discussion. Chaotic 
inflationary models, associated with a large energy (~ GUT scale) of y^/^ ~ lO^^GeV, 
predict a very large number of inflationary e-folds, ^ 62. Long-lasting inflation erases 
all classical anisotropics and inhomogeneities of the pre-inflationary stage. However, scalar 
and tensor vacuum fluctuations during inflation lead to almost scale free post-inflationary 
scalar and tensor metric inhomogeneities around our smooth observable FRW patch. 

In particular, the amplitude of the gravitational waves generated from the vacuum fluc- 
tuations during inflation is proportional to V^^"^, — ^^Hj ~ (where Mp is the reduced 
Planck mass). There are signiflcant efforts to measure the i?-mode of AT/T polarizations, 
since this will provide a direct probe of the scale of inflation. The current 95% CL. limits on 
r (ratio of the tensor to scalar amplitudes of cosmological fluctuations) r <, 0.43 (WMAP- 
only) and r 0.2 (WMAP plus acoustic baryon oscillation, plus supernovae) |l| shall be 
improved to r <, 0.1 by the Planck mission j^, to r <, 0.05 by the Clover 3], EBEX {4], 
and Spider [^J experiments (see 16] for the study of a mission that can improve over these 
limits). While these limits imply a detection in the case of high energy inflation, a number of 
other inflationary models, including many of the string theory constructions have lower en- 
ergy, and therefore lead to GW of much smaller amplitude, which are virtually unobservable 
through B mode polarization ^. 

In anticipation of the null signal observation of the primordial GW from inflation, it 
is worth thinking about other implementations of this result for the theory of inflation, 
besides putting limits on the energy scale V^^'^. There are models of inflation (including 
many string theory inflationary models) where the total number of e-folds, A^, does not 



^ Future gravitational waves astronomy may allow to probe r up to the level 10 ^ with BBO [?[ or ultimate 
DECIGO 8} direct detection experiments. 
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exceed the minimum ([I]) by a large number. If the extra number of e-folds AA^ beyond 
([1]) is relatively small then pre-infiationary inhomogeneities of the geometry are not erased 
completely, and their residuals can be subject to observational constraints. In the context of 
this idea, in this paper we suggest an additional mechanism to have observable gravitational 
waves associated with inflation. These gravitational waves are very different from the GW 
generated from the vacuum fluctuations during inflation. Firstly, they are the residual tensor 
inhomogeneities from the pre-inflationary stage. Secondly, they can be of a classical, rather 
than quantum, origin. Thirdly, while their initial amplitude and spectrum are given by 
the initial conditions, they are significantly affected by the number of "extra" e-folds AA^. 
Therefore, observational limits on gravity waves result in constraints on a combination of 
AA^ and of the initial amplitude. 

The choice of the initial geometry of the universe before inflation is wide open. In 
principle, one may assume an arbitrary geometry with signiflcant tensor inhomogeneities 
component, and much smaller scalar inhomogeneities. This choice is, however, very artiflcial. 
A much more comfortable choice of the pre-inflationary stage will be a generic anisotropic 
Kasner-like geometry with small inhomogeneities around it. The origin of the anisotropic 
universe with the scalar fleld can be treated with quantum cosmology, or can be embedded 
in the modern context of the tunneling in the string theory landscape. In fact, a Kasner-like 
(Bianchi I) space was a rather typical choice in previous papers on pre-inflationary geometry, 
see e.g. [9|]. Most of the works on an anisotropic pre-inflationary stage aimed to investigate 
low the initial anisotropy is diluted by the dynamics of the scalar fleld towards inflation 

The formalism of linear fluctuations about an anisotropic geometry driven by a scalar 
fleld toward inflation was constructed only recently H, 12, Q, [l^. Besides the technical 



aspects of calculations of cosmological fluctuations, there is a substantial conceptual differ- 
ence between computations in the standard inflationary setting and in the anisotropic case. 
For an isotropic space undergoing inflationary expansion, all the modes have an oscillator- 
like time-dependence at sufficiently early times, when their frequency coincides with their 
momentum. One can therefore use quantum initial conditions for these modes. This is no 
longer the case for an expansion starting from an initial Kasner singularity. In this case, a 
range of modes, which can potentially be observed today (if AA^ is not too large), are not 
oscillating initially and therefore cannot be quantized on the initial time hyper-surface; as 
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a consequence, there is an issue in providing the initial conditions for such modes. For this 
reason we will adopt another perspective, namely, we will consider generic small classical 
inhomogeneities around the homogeneous background, as an approximation to the more 
generic anisotropic and inhomogeneous cosmological solution. 

Equipped with this philosophy, we consider an anisotropic expanding universe filled up by 
the scalar field with a potential V{(j)) which is typical for the string theory infiation. We add 
generic linear metric fiuctuations about this geometry. The evolution of these fiuctuations 
is by itself an interesting academic subject. However, it acquires a special significance 
in the context of the GW signals from infiation, because of a new effect that we report 
here of amplification of long-wavelength GW modes during the Kasner expansion. This 
growth terminates when a mode enters the "average" Hubble radius (the average of that 
for all the three spatial directions), or, for larger wavelength modes, when the background 
geometry changes from anisotropic Kasner to isotropic infiationary expansion. We perform 
explicit computations in the case of an isotropy of two spatial directions. In this case the 
computation becomes much more transparent and expli citly k dependent. Fluctuations for 



arbitrary a, b, c were considered in the formalism of 



where the k dependence is not 

explicit. We verified that our results agree with |12l . |14| | in the axisymmetric b = c limit. 
We find that only one of the two GW polarizations undergoes significant amplification. 
Therefore, even if we assume for simplicity equi-partition of the amplitudes of the three 
inhomogeneous physical modes of the system (the scalar and the two GW polarization) at 
the initial time, the final spectra that will be frozen at large scales in the infiationary regime 
will be very different from each other, in strong contrast to what is obtained in the standard 
inflationary computations. 

This result can have different consequences, that we explore in the present work. Suppose 
that the growing GW mode is still linear (but signiflcantly exceeds other modes) when the 
space becomes isotropic. Then, we can have signiflcant yet linear classical GW fluctuations 
at the beginning of inflation, say of amplitude ^ C>{1). If the modes which correspond to 
the largest scales that we can presently observe left the horizon after the flrst AA^ e-folds of 
inflation, their amplitude decreased by the factor e~^^ in this period. If AA^ is relatively 
small, say ~ C(10) the freeze out amplitude of these GW modes would be ~ O(10~^). 
The angular spectrum of these GW will rapidly decrease as the multipole number i grows, 
since smaller angular scales are affected by modes which spend more time inside the horizon 
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during the inflationary stage. 

Suppose instead that the growing GW mode becomes non-hnear before the onset of 
inflation. In this case the background geometry departs from the original onset. 

Besides the phenomenological signatures, it is interesting to study the origin of the am- 
plification of the GW mode. It turns out that the effect of GW amplification is related 
to the anisotropic Kasner stage of expansion. Therefore we will separately study GW in 
the pure expanding Kasner cosmology. For completeness, we also include the study of GW 
in a contracting Kasner universe, which is especially interesting due to the universality of 
anisotropic Kasner approach to singularity. 

The plan of the paper is the following. In Section [TTl we discuss the evolution of the 
anisotropic universe driven by the scalar field towards inflation. In Section IIIII we briefly 
review the formalism of the linear fluctuations in the case of a scalar field in an anisotropic 
geometry, paying particular attention to the GW modes. In Section [IV] we compute the 
amplification of one of the two GW modes that takes place at large scales in the anisotropic 
era. In Section |V] we discuss instead the evolution of the other two physical modes of 
the system. In Section IVII we study the evolution of the perturbations in a pure Kasner 
expanding or contracting Universe. In Section IVlII we return to the cosmological set-up, and 
we compute the contribution of the GW polarization amplified during the anisotropic stage 
to the CMB temperature anisotropics. In particular, by requiring that the power in the 
quadrupole does not exceed the observed one, we set some limits on the initial amplitude of 
the perturbations vs. the duration of the inflationary stage. In Section rVIIII we summarize 
the results and list some open questions following from the present study, which we plan to 
address in a future work. 

II. BACKGROUND GEOMETRY 

The anisotropic Bianchi-I geometry is described by 

ds"^ = -dt^ + a^dx^ + l?dy'^ + c^dz^ , (2) 

where a(t), 6(t), c{t) are the scale factors for each of the three spatial directions. 

We consider a scalar field in this geometry. Many string theory inflationary models (for 



examples see 



13, 



18j) have a very flat inflationary potential which changes abruptly around 
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its minima. Therefore, to mimic this situation, we will use a simple inflaton potential 



(3) 



which has quadratic form y ~ 0^ around the minimum, and is almost flat V ^ Vq away 
from it. To obtain the correct amplitude of scalar metric perturbations from inflation, we 
set 



^ = 10-^ 



Vr 



1/4 



lO^^GeV . 



(4) 



The background dynamics is governed by the Einstein equations for the scale factors in 
the presence of the effective cosmo logical constant Vq, plus a possible contribution from the 
kinetic energy of the scalar field. Quantum cosmology or tunneling models of the initial 
expansion favor a small scalar field velocity. Therefore we select the small velocity ini- 
tial conditions 0^ <^ Vq. In this case the generic solutions of the Einstein equations with 
cosmological constant for a, b, c are known analytically (see e.g. 1J|) and can be cast as 



a(t),b{t),c{t)] = (a 



I bin 1 Cj' 



^171"! ^171"! ^in ) 



sinh(3i7ot) 



1/3 



tanh I 2-^0^ 



nK-i/3 



(5) 



where are the Kasner indices, X]i=iPi — X]i=iPi — ^i ^^'^ -^o is the characteristic 
time-scale of isotropization by the cosmological constant, = while 
the normalizations of the three scale factors 

For earlier times t <^ I/Hq the anisotropic regime is described by the vacuum Kasner 
solution 

(a(t),6(t),c(t)) = {ao,bo,Co)-tP^ , (6) 

1 /3 

which corresponds to an overall expansion of the universe (the average scale factor (a6c) = 
t^/^ is increasing), although only two directions are expanding (two positive Pi-s) while the 
third one is contracting (the remaining pi is negative). 

For later time t ^ 1/Hq the universe is isotropic and expanding exponentially 



a(t),6(t),c(t) ) = (ao, 6o,co) ■ e 



Hot 



(7) 

where the constant normalizations ao, &05 Co are typically chosen to be equal. 

It is instructive to follow the evolution of the curvature in the model. The Ricci tensor 
is (almost) constant throughout the evolution up to the end of inflation 

1 



R = 12m 



(8) 



At earlier times the Weyl tensor - i.e. the anisotropic component of the curvature tensor - 
gives 

C^.C'-C,„^ = -15i!l|«5, (9) 

and, for t <^ I/Hq, \fC^ is much bigger than the isotropic components ([8]). This is why ini- 
tially the contribution from the effective cosmological constant is negligible, and the vacuum 
Kasner solution ([6]) is a good approximation. In contrast, at later times l/-f^o 

C^^'^'C^,,, ~ e-^^°* , (10) 

and the anisotropic part of the curvature becomes exponentially subdominant relative to its 
isotropic part driven by the cosmological constant. This is an illustration of the isotropiza- 
tion of the cosmological expansion produced by the scalar field potential. The timescale for 
the isotropization is tiso = l/-^o • 

In the following Sections we will study the equations for the linear fluctuations around 
the background (EI), ([5]). These equations become significantly simpler and more transparent 
for the particular choice of an axi-symmetric geometry e.g. when c = b, and the metric is 

ds^ = -dt^ + a^dx^ + b^ {dy^ + dz^) . (11) 

While the effect we will discuss is generic, for simplicity we will adopt the simpler geometry 
f fTTj) rather than the general Bianchi-I space ([2]). In this case, the early time solution is a 
Kasner background with indices pi = —1/3,^2 = Ps = 2/3 ^. Also, it will be useful to define 
an "average" Hubble parameter H and difference h between the expansion rates in x and y 
(or z) directions as 

jj^lh±2lh^ h^'ksJk^ ff. = ^ff, = ^ (12) 

3 V3 a b 

At earlier times Ha ^ Hh ^ while, at late times. Ha = Hf, ^ Hq. The equation for 
the homogeneous scalar field is 

+ 3H^ +V^ = . (13) 



^ There is another axisymmetric asymptotic Kasner solution, with indices pi — 1, P2 = Ps = 0- The solution 
is a very special one, since it is the only Bianchi-I model with cosmological constant that is regular at 
t = Q (as can be easily checked by computing the curvature invariants; e.g. = ISTJg at i = 0). For 
Hq — 0, this space is actually Minkowski space-time in an accelerated frame. Due to its special nature, 
we disregard this solution in the present study. 
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Since the value of H is very large initially, if ~ ^, the Hubble friction keeps the field 
(practically) frozen at = 0j„ during the anisotropic stage. For t > t-^^o, the universe 
becomes isotropic, and it enters a stage of slow roll inflation until (j) rolls to the minimum 
of its potential. 

We also will use another form of the metric (|TT1) . with the conformal time rj. There is 
ambiguity in the choice of 77, related to possible different choices of the scale factors in its 
relation with the physical time. We will use the average scale factor 

a,^^{ah''Y'\ (14) 

and define t) through 

dt= [ab^y^'^dr] , (15) 

which, at early times, gives i] oc t^^^. In this variable, the line element (fTTj) reads 

2/3 



ds^ = {ab^) 



2\2/3 



-d?]"^ + (^^^ dx^ + [dy'^ + dz^) 



(16) 



In the following, dot denotes derivative wrt. physical time t, and prime denotes derivative 
wrt conformal time. Moreover, we always denote by Ha and the Hubble parameters with 
respect to physical time. 

III. LINEAR FLUCTUATIONS 

In the FRW universe with a scalar field there are three physical modes of linear fluctua- 
tions. Two of them are related to the two polarizations hj^ and of the gravitational waves, 
and one to the scalar curvature fluctuations v induced by the fluctuations of the scalar field 
5(f). All three modes in the isotropic case are decoupled from each other. The formalism for 
the linear fluctuations on a FRW background has been extended to the Bianchi-I anisotropic 
geometry in 
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id, [isl, [3]. Again, there are three physical modes; however, in the general 



case of arbitrary a, 6, c the modes are mixed, i.e. their effective frequencies in the bi-linear 
action are not diagonal, as it is the case in the isotropic limit. 

In a special case of the axi-symmetric Bianchi I geometry flTT]) one of the three linear 
modes of fluctuations, namely, one of the gravity waves modes, is decoupled from the other 
two. This makes the analysis of fluctuations much more transparent than the general a 7^ 
b ^ c case. While the effects we will discuss, we believe, is common for arbitrary a, 6, c, 
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we will consider linear fluctuations around the geometry (fTT]) . The computation follows the 



formalism of 



where the reader is referred for more details. 



The most general metric perturbations around fill I) can be written as 



(17) 



f -al, (1 + 2 $) a diX a,, b (5 , + Bi) ^ 

a2(l-2^) abdi (^B^i + R^^ 
\ 62[(l_2S)5,^. + 2E,,+%,)] J 

where the indices i,j = 2,3 span the {y, z} coordinates of the isotropic 2d subspace. The 
above modes are divided into 2d scalars ($, x, B, \E', B, S, E) and 2d vectors {Bi, Bi, Ei, 
subject to the transversality conditions diBi = diBi = diEi = 0) according to how these 
modes transform under rotations in the isotropic subspace. The two sets of modes are 
decoupled from each other at the linearized level. In addition, there is the perturbation of 
the inflaton field = (j) + 5(j) , which is also a 2d scalar. 
The gauge choice 

5gii,2ds = 5gij = ^ , (18) 

corresponding to B = 'L = E = Ei = {], completely fixes the freedom of coordinate 
reparametrizations. It is convenient to work with the Fourier decomposition of the linearized 
perturbations. We can therefore fix a comoving momentum k, and study the evolutions of 
the modes having that momentum. Since modes with different momenta are not coupled at 
the linearized level, this computation is exhaustive as long as we can solve the problem for 
any arbitrary value of k. More precisely, we denote by the component of the momentum 
along the anisotropic x direction, and by kx the component in the orthogonal y — z plane. We 
denote by and by px the corresponding components of the physical momentum. Finally, 
we denote by k and p the magnitudes of the comoving and physical momenta, respectively. 
Therefore, we have 

k' = kl + kl , P'=Pi+PT=(^)'+(y)' (19) 

To identify the physical modes, one has to compute the action of the system up to 
the second order in these linear perturbations. One finds that the modes ^, x, B, and Bi, 

^ Notice that transverse 2d vectors have one degree of freedom; contrary to the 3d case, there are no 
transverse and traceless 2d tensors. 
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corresponding to the Sgofj, metric perturbations, are nondynamical, and can be integrated out 
of the action. This amounts in expressing the nondynamical fields (through their equations 
of motion) in terms of the dynamical ones, and in inserting these expressions back into the 
quadratic action. For instance, for the nondynamical 2d vector mode one finds 



!,\ 1/3 2 / 

b\' pi ^ 



(20) 



The analogous expressions for the 2d nondynamical scalar modes can be found in ISj. 

In this way, one obtains an action in terms of the three remaining dynamical modes 
6(j), ^E', -Bj. Once canonically normalized, these modes correspond to the three physical per- 
turbations of the system. The canonical variables are 



V 



+ 



Pt' 



HaPl + H,{2pl+pl) 



a/2 aavMpp^iffc 
HaP% + H,{2pl+p%) 



and 



-TZ X — /-- ^scv ^ij Pi -Dj 

V2 P 



(21) 



(22) 



where eij is anti-symmetric and ei2 = 1 (we stress that encodes only one degree of 
freedom, since, due to the transversality condition of the 2d vector modes, PiBi = 0). 

The dynamical equations for the modes if+ and V are coupled to each other, while that 
of the Hx mode is decoupled 








LOT, CJT 




'11 ^12 
a2 '^22 

The explicit expressions for the frequency matrix uj, 
given in jiaij . 

In the limit of isotropic background, b ^ a, also the 2d scalars decouple, and the fre- 



- . (23) 
1,2) are rather tedious and 



quencies become 



X' ^^22 



11 



12 



(24) 



Therefore, the mode V becomes the standard scalar mode variable 



19| , associated to the 



curvature perturbation, while the modes if+ and i^x are associated to the two polarizations 
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of the gravitational waves. Also notice that without the scalar field there are two physical 
modes which, due to the residual 2d isotropy, are decoupled. 



IV. DECOUPLED TENSOR POLARIZATION 

Since there are no vector sources, the 2d vector system describes a polarization of a 
gravitational wave obeying the "free field equations" 5Rfj,u = (which reproduce the equation 
fl20|) and the first of fl23|) ). As we now show, this mode undergoes an amplification, which 
does not occur for the other two modes. This can be understood from considering the 
frequency of this mode 



^^-Pl+Pt + ^ + ^-^^'^-^^) • (25) 

,2 



In the isotropic case, for which uj\ is given by fl2^ . each mode is deeply inside the horizon, 
k/a ^ H, at asymptotically early times. This is due to the fact that H is nearly constant, 
while k/ais exponentially large at early times. As a consequence, — in the asymptotic 
past, and the mode oscillates with constant amplitude. In the present case instead 

PL oc a-^ oc r^^/^ ^ t^/^ , pt oc 6"^ oc r/"^ ^ t'^/^ , Ha, oc r^'^/^ ^ ^-i ^26) 

at early times (r^, t —>■ 0^). Namely, as we go backwards in time towards the initial sin- 
gularity, the anisotropic direction becomes large, and the corresponding component of the 
momentum of the mode is redshifted to negligible values. On the contrary, the two isotropic 
directions become small, and the corresponding component of the momentum is blueshifted. 
However, the magnitude of the two Hubble parameters increases even faster. Therefore, 
provided we can go sufficiently close to the singularity, the early behavior of each mode is 
controlled by the negative term proportional to the Hubble parameters in eq. (!25|) . 

To be more precise, if we denote by ag and bo the values of the two scale factors at some 
reference time rjo close to the singularity, we have 



(27) 



flav (^o) T 

where the first term in the expansion comes from the terms proportional to ^ in eq. 
while the second term from the component of the momentum in the isotropic plane (cf. the 
early time dependences with those given in eq. (I26l) ). We see that the frequency squared 
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is negative close to the singularity, so that the mode H-^ experiences a growth ^. In a 
pure Kasner geometry, the relations fl26|) hold at all times. Therefore, one would find 
Pl ^ Pt ^ Ha^h at asymptotically late times. For brevity, we will loosely say that the 
mode "enters the two horizons H~l" at late times; the meaning of this is simply that the 
frequency is controlled by the momentum in this regime, ~ a^-^p"^ > 0, and the mode 
Hx enters in an oscillatory regime. 

This simple description is affected by two considerations: firstly, we do not set the initial 
conditions for the modes arbitrarily close to the singularity, but at some fixed initial time 
TjQ ; Secondly, the geometry changes from (nearly) Kasner to (nearly) de Sitter due to the 
infiaton potential energy. Consequently, there are three types of modes of cosmological size. 
I: Modes with large momenta start inside the two horizons at r/o • They oscillate (tu^ > 0) all 
throughout the anisotropic regime, and they exit the horizon during the infiationary stage. 
II: Modes with intermediate momenta, for which (!27|) is a good approximation at tiq. These 
modes enter the horizons, and start oscillating, at some time rj > t]o during the Kasner era; 
they exit the horizon later during infiation. Ill: Modes with small momenta, that are always 
outside the horizons, and never oscillate during the Kasner and infiationary regimes. 

These considerations are crucial for the quantization of these modes. We can perform 
the quantization only as long as cu^ > 0, and the mode is in the oscillatory regime. As we 
mentioned, during infiation, this is always the case in the past. Moreover, the frequency 
is adiabatically evolving {u' ^ cj^), and one can set quantum initial conditions for the 
mode in the adiabatic vacuum. This procedure is at the base of the theory of cosmological 
perturbations, and results in a nearly scale invariant spectrum at late times, once the modes 
have exited the horizon and become classical. In the case at hand, we cannot perform this 
procedure for modes of small momenta / large wavelength (modes III above) . If infiation lasts 
sufficiently long, such modes are infiated to scales beyond the ones we can presently observe, 
so that the inability of providing initial quantum conditions is irrelevant for phenomenology. 
However, if infiation had a minimal duration, this problem potentially concerns the modes 
at the largest observable scales. 

Irrespectively of the value of the frequency, it is natural to expect that the modes possess 

^ A tachyonic frequency does not necessarily imply a growth of the physical fluctuations. For instance, this 
does not happen for the GW modes in the FRW geometry. 
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some "classical" initial value at t]o . In the following we discuss the evolution of the pertur- 
bations starting with these initial conditions. Although we do not have a predictive way to 
set these initial values, we can at least attempt to constrain them from observations. In ad- 
dition, we should worry whether the growth of Hy. can result in a departure from the Kasner 
regime beyond the perturbative level, in which case the background solution described in 
the previous Section may become invalid (we discuss this in Section rvTj) . 

As long as the frequency is accurately approximated by fl27|) in the Kasner regime, the 
evolution eq. for the mode Hy. is approximately solved by 

1/3 



H^=Ciy/ri J3 



> 1/3 
ao ^ 



2 ( ^ ) 



(28) 



where Ci,2 are two integration constants, while J and Y are the Bessel and Neumann func- 
tions. The first mode increases at early times (small argument in the Bessel function) oc 1]'^, 
while the second one decreases as t]~^ . We disregard the decreasing mode in the following 
computations, C2 = (moreover, this mode diverges at the singularity). 

Rather than the time evolution of , we show in Figure [1] that of the corresponding 
power. We note the very different behavior obtained for large (III), intermediate (II), and 
small (I) scale modes. We choose to define the power spectrum as 



2 

p 



p'\H^r (29) 



This definition coincides with the standard one (see for instance [20|) as the universe be- 
comes isotropic. In particular, the power spectrum is frozen at large scales in the isotropic 
inflationary regime. Clearly, there is an ambiguity in this definition at early times, when 
the two scale factors differ (there is no a-priori reason for the choice of the average scale 
factor in this definition). This arbitrariness affects the time behavior shown in the Figure. 
However, it does not affect the relative behavior of the large vs. intermediate vs. small scale 
modes. Moreover, if we analogously define the power spectra for the 2d scalar modes, the 
relative behavior of these two types of mode (Figure [1] vs. Figure H] ) is also unaffected by 
this arbitrary normalization. 

The reality of this instability is demonstrated in Section Wl\ where we compute the 
squared Weyl invariant due to these fluctuations (more precisely, we do so in an exact Kasner 
background, which, as we have remarked, coincides with the cosmological background at 
asymptotically early times). 

13 



1e+15 
1e+10 
100000 
1 

of 1e-05 
of 1e-10 
1e-15 
1e-20 
1e-25 
1e-30 

1e-08 1e-06 0.0001 0.01 1 

t/tiso 

FIG. 1: Power in different modes of the decoupled tensor polarization, normalized to its initial value 
{ki = kx for all the cases shown). We show the time evolution starting from /iq = —1^^\/Vq/M 
(defined in eq. [T2]) . The universe isotropizes at t\so = 1/^iso- The quantity /ciso is defined to be 
the momentum of the modes exiting the horizon at this time. Large scale modes do not exhibit 
oscillatory behavior, and grow during the anisotropic stage. Intermediate scale modes enter the 
horizons during the anisotropic era; this terminates their growth. Small scale modes are in the 
oscillatory regime all throughout the anisotropic phase, and do not experience any growth. All 
modes shown freeze during the inflationary stage. Notice the very different final values obtained 
in the three cases shown. 

In Figure [2] we show the power spectrum (normalized to the initial value for each mode) 
for the same background evolution as in the previous Figure, at some late time during 
inflation, when all the modes shown are frozen outside the horizon. As the approximate 
solution (!28|) indicates, the growth of iJx occurs as long as the transverse momentum px is 
smaller than the Hubble rates \Ha,b\ ■ Therefore, modes with smaller kx experience a larger 
growth. We denote by 6 the angle between the comoving momentum, and the anisotropic 
direction, 

kL = kcos9 , kT = ksm9 (30) 
14 




^so 



Therefore, in general, we expect a greater growth at smaller values of k (for any fixed 6) and 
at smaller values of 6 (for any fixed k). ^ This behavior is manifest in Figure [2l Modes with 
k <^ fciso experience the same growth during the Kasner era (since the leading expression 
for uj"^ is independent of the momentum in this regime). Then, the modes shown in the 
Figure exhibit a very strong ~ k~^ dependence for ki^o ^ k < few x 10^ ki^o ■ We also see 
an increase of the power as 6 decreases. We stress that Figure [2] shows the contribution of 
each mode to the power spectrum normalized to the value that that contribution had at the 
initial time r]Q . Therefore, if the original spectrum has a scale, or an angular dependence, 
this will modify the final spectrum (for comparison, for the isotropic computation of modes 
with adiabatic quantum initial condition, Ph^ oc k^ |-f^xP oc k'^ at early times). 

1e+10 fznn^nninnii^. — • • — ^ • • — ^-i • • — ^ • • — n 
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1e-25 I — ■ ■ — ^ — ■ — ■ ■ — ^ — ■ — ■ ■ — ^ 

0.1 1 10 100 1000 10000 

k/k|5o 

FIG. 2: Amplification of the power spectrum between the initial time (same as in the previous 
Figure) and some time during inflation for which the modes shown are outside the horizon, and 
the value of Ph^ is frozen. As explained in the main text, the growth is greater at small values of 
k and 9. 

The large growth at small 6 is more manifest in Figure [3l The smallest angles shown 
in the Figure correspond to kx <^ ki, but to ^ Pl at the initial time (this is due to 

^ Notice that, due to the planar symmetry of the background, the same results are obtained at and n — d. 
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the different behavior of the two scale factors in the anisotropic era). In this region, the 
spectrum exhibits a milder Ph^, oc dependence than for intermediate values of 6. Finally, 
one may also consider smaller angles than those shown in the Figure, for which < Pl 
initially. We have found that final spectrum becomes 6'— independent in this region. 

1e+15 
1e+10 
100000 
1 

o 
Q. 

~~^x 1 e-05 
1e-10 
1e-15 
1e-20 
1e-25 

0.0001 0.001 0.01 0.1 1 

sin 

FIG. 3: Angular dependence of the power of different modes. We show the value reach by the 
mode once it is frozen in the inflationary stage, normalized to the initial value (the background 
evolution is identical to that of the previous two Figures) . In the region 9 <^ 1 shown in the Figure 
{kj- <C k^) the final value exhibit a 1/ sin 6 dependence. 

We conclude this section by discussing how the growth scales with the initial time. As 
long as \Ha,b\ ^ Pt ^ Pl, the power ( 129|) of a mode grows as r]^^'^ (as can be easily seen 
by combining the time dependences oc r^^/^, p ~ oc ?7, Hy^ (x rf). We use the initial 
value of h (the difference between the two expansion rates, defined in eq. (fT2|) ) as a measure 
on the initial time, since, contrary to the conformal time, this quantity is not affected by 
the normalization of the scale factors. Starting with a greater value of |/io| corresponds 
to starting closer to the initial singularity, and, therefore, to a longer phase in which Ph^ 
grows. Since h oc t^^ oc r/^'^/^ in the Kasner regime, the ratio Ph^/Pq oc \hQ\ in the region 
in which the growth takes place. Although we do not show this here, we have verified that 
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this scaling is very well reproduced by the numerical results. 



V. COUPLED TENSOR POLARIZATION-SCALAR MODE PAIR 



As discussed in Section [nil the two other physical modes of the system are coupled to each 
other in the anisotropic era. The evolution equations for the coupled system are formally 
given in fl23l) . At early times, we find 



1 r/0 _^ 



^12 



^21 



(31) 



Therefore the coupling between the two modes can be neglected also at asymptotically 
early times. The main difference with the analogous expression for the decoupled tensor 
polarization, eq. ([27]), is that the squared eigenfrequencies of the two modes are now positive 
close to the singularity; therefore the two modes V, do not experience the same growth 
as Hx . Indeed, as long as the expressions ( |3T1) are good approximations, we find the solution 



1/3 



bo 



1/3 



VVoV 



(32) 



and an identical one for V with the replacement ^ C*! 2 integration constants. 

Close to the singularity, the two modes grow as y/rj and ^Ini] , respectively. Analogously 
to what we did for the Hx polarization, we disregard the mode which grows less at early 
times, Cf + = C^ = 0. 

We define the power spectra for the tensor polarization, and for the comoving curvature 
perturbation R with the same prefactor as Ph^, cf. eq. f l29|) . 



"72" 



R 



H 



''av 3 



(33) 



TT- V y 27r2 

We see that, contrary to what happened for , the coupled perturbations, and the cor- 
responding power spectra, do not grow while outside the horizons during the anisotropic 
era. 

This effect is also manifest in Figure [21 where we show the spectra of the tensor mode 
iJ+ and the comoving curvature for the same range of momenta, and for the same angles 6, 
as those of Ph^ shown in Figure [2l 



17 



1e-05 i 



1e-10 - 



1e-15 



1e-20 



1e-25 



1e-30 




k = 0.01 k:„ 
k= 10^ k:„ 



g. 



1e-08 1e-06 0.0001 0.01 



1 

1e-05 
1e-10 
1e-15 
1e-20 
1e-25 
1e-30 
1e-35 
1e-40 



1e-08 




k = 0.01 k: 
k = 1 k| 
k = 1 0* k; 



1e-06 



FIG. 4: Contribution of different modes to the power spectrum of the tensor polarization (left 
panel) and of the comoving curvature perturbation R (right panel), normalized to its initial value. 
The background evolution, and the momenta shown, are the same as in Figure [TJ Contrary to the 
decoupled tensor mode, shown in Figure [H the spectra do not grow while the modes are outside 
the horizons in the anisotropic regime. 

VI. INSTABILITY OF KASNER SOLUTION AGAINST GRAVITATIONAL 
WAVES 



The main result of the previous section was a significant amplification of the mode H^., 
compared to the milder amplification of the mode in the anisotropic background which 
is undergoing isotropization due to the effect of a scalar field. This growth can be ascribed to 
the instability of the Kasner geometry, either contracting or expanding, against gravitational 
waves which we are going to report in this Section. Therefore in this Section we consider 
linearized gravity waves around an expanding and a contracting Kasner solution, without 
the presence of the scalar field, nor its fluctuation. In this case there are only the two 
decoupled modes and if+ . 

The claim of instability of the Kasner solution against the growth of the GW sounds at 
first glance heretic, at least for the contracting branch, in the light of the universality of the 
Belinskii-Khalatnikov-Lifshitz (BKL) oscillatory regime of the Kasner epochs approaching 
the singularity. In fact, it is not, and, on the contrary, it is compatible with the BKL 
analysis. Moreover, our finding of the GW instability suggests a new interpretation of the 
phenomena connected to the instability, discussed by BKL and others for the contracting 
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FIG. 5: Decrease of the power spectra of the tensor polarization (left panel) and of the comoving 
curvature R (right panel). The Figure shows the ratio between the power spectra at some time 
during inflation (when the modes are outside the horizon in the isotropic regime, and the spectra 
are frozen), and the initial time. The decrease should be compared with the growth experienced 
at large scales by Pf/ , shown in Figure [2j 
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35]. 



Kasner geometry in very different formalism and language 

In this Section we first perform linearized calculations for the classical gravitational waves 
around expanding and contracting Kasner solutions, and demonstrate their instabilities in 
terms of the evolution of the Weyl tensor invariant C^^'^'^Cf^ppa-, which is independent of the 
gauge choice. We then connect the result with the BKL analysis. 

The background line element is given by equation (fT6|) . with the scale factors 

1/2 



a = ao 



V 



b = bo 



(34) 



This compact notation describes two disconnected geometries, at negative and positive con- 
formal times, respectively. The algebraic expressions below simplify if we introduce the time 
rj^ in which the normalization of the two scale factors coincide. Therefore, rather than (!34ll . 
we can also use 



a = a* 



V 



1/2 



The two Hubble rates are 



2 a* 



1/2 1 



1/2 1 
7] 



(35) 



(36) 
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while the average scale factor is = 
related by 

dt = ttav dr] 



■ The physical and conformal times are 

(37) 



t 



2 a* 



V* 



1/2 



V 



For future use, we also define to be the physical time corresponding to 77=,, . 

The solution with negative conformal time describes an overall contracting geometry, 

1 /2 

c^av oc {—Tj) , which crunches into the singularity at 77 = 0. The solution with positive 
conformal time describes instead an overall expanding space, a^v oc 77^^^, originating at the 
singularity at rj = . As in the previous Sections, we restrict the computation to the simpler 
case of a residual 2d isotropy between two spatial directions. We expect that the instability 
occurs for general Kasner indices (pi, p2, Ps). 

Now we turn to the linearized perturbations satisfying the vacuum equations SRf^i, = 0. 
There are two gravitational waves polarization perturbations. We consider a single mode 
with a given comoving momentum with components ki and /c^. We denote by ^ the vector 
in the y — z plane along the direction of kx, 



H (r/, x) = Q-^k^^-^'^T-iH (77, fci, kr) + h.c. 
The two GW polarizations obey the equations 



m 







H'i + ul H. 







(38) 



(39) 



where the effective frequencies can be written in compact form 



2pI + 20pIp^t-8p^t 



P' + Hi 



p" 



,16pI + 296pIpI+P^t 



(40) 



{Apl + p%f 

both for the contracting and the expanding backgrounds. We recall that the physical mo- 
menta are related to the comoving one by the relations (fT9l) . 

After solving the two equations fl39|) . we can compute the metric perturbations f|T7j) and 
the Weyl tensor of the background plus perturbations. The square of the Weyl tensor, once 
expanded perturbatively, has the following schematic structure 



flupcr 



+ C5C + + 



(41) 



where is the Weyl invariant for the non-perturbed background solution, C SC is the 
term linear in and H^, and 60"^ is the term quadratic in the perturbations. We do the 
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computation for the two polarizations separately. For instance, for the 2d vector modes, we 
compute the Weyl tensor in terms of the background and of the metric perturbations B3 
and B3. We then relate the two perturbations to Hy. through Eqs. (HUl) and expressing 
spatial derivatives in terms of the comoving momenta, see Eq. (1551) . In this way, we can 
write the expression (HT]) in terms of H^. , if* , and their time derivatives. Finally, we insert 
in this expression the solutions of Eq. (139|) . The procedure for the mode if+ is analogous. 

We compare the second and third term in fj4T]) with the background term. A growth of 
the ratios CdC/C"^ (denoted as = 5C/C)^ or dC^/C"^, signals an instability of the Kasner 
geometry. The background (zeroth order) Weyl invariant is 

^ - atv'~ 27 • ^^^^ 
The first order term C 5C vanishes identically for the polarization. It is nonzero for H^, 



and it oscillates in space as exp 



The second order term for the 



± i {kL X + kT ■ tj 

Fourier modes or has a part which is constant in space, plus two parts that oscillate 
in space as exp ± 2i {ki x + kx ■ - In the following, we disregard the oscillatory parts 
in (5C2 . 

The solutions and are either monotonically evolving in time, or oscillating, with 
an envelope that is monotonically evolving in time. The oscillatory regime takes place when 
the mode has a wavelength shorter than the Hubble radii, p > \Ha\, \Hb\, and are absent 
in the opposite regime. Eqs. fl26l) show the time dependence of the momentum and of the 
Hubble rates for the expanding Kasner solution. We see that, if we consider a complete 
background evolution, any mode starts in the large wavelength regime, and then goes in the 
short wavelength regime. Therefore, we expect that a mode is in the non-oscillatory regime 
sufficiently close to the singularity, and in the oscillatory regime sufficiently far from it (this 
behavior is manifest in the time evolutions shown). The same is true for the contracting 
background solution (with the obvious difference that early and late times interchange). 

The time dependence of the Weyl invariant (or of its amplitude, when the mode is os- 
cillating) is summarized in Table 1. More specifically, we show the ratio between the terms 
proportional to the perturbations and the background one, both for the expanding and the 
contracting Kasner. In the expanding case, a mode evolves from the large scale to the short 
scale regime. The opposite happens in the contracting case. Notice the the short scale 
behaviors in the expanding and contracting cases coincide. The same is not true for the two 
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EXPANDING 

(0 < t < oo) 



CONTRACTING 

(-00 < t < 0) 



5C/C for 





dC^/C^ for 





bC^jC^ for 



large scales 
short scales 



short scales 
large scales 



TABLE I: Summary table for the time dependence of 6C /C and of the (non-oscillating in space ) 
part of (5C72/C72. 

large scales behaviors. The reason is that, in the expanding case, we set to zero one of the 
two solutions of Eq. ( !39l) that is decreasing at early times, and that would diverge at the 
initial singularity. The different behaviors are discussed in details in the next Subsection 
and in Appendix \M 



A. Contribution of the i^x mode to the Weyl invariant 

In this Subsection, we compute the contribution of the Hy^ to the square of the Weyl ten- 
sor, both for an expanding and a contracting Kasner geometry. The analogous computation 
for the mode if+ can instead be found in Appendix [Al 



1. Hx mode during expansion 

Plugging fl35l) into fj40|) . the frequency has the large scales (early times) and short scales 
(late times) expansions 



1 + 'L* 



large scales 
short scales 



(43) 



Consequently, we have the asymptotic solutions 



Ci y/v J3 (2 y/T]^ rj) , large scales 



Ql. e* 2j7, _j_ g ' 2j7: 



C2 



short scales 



(44) 



where Ci, Ci, C2 are three integration constants. In the early time solution we have dis- 
regarded a decaying mode that would diverge at 77 — , and where the expression in the 
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second line is the large argument asymptotic expansion of the parabolic cylinder functions 




(45) 



which are solutions of the evolution equation with the short scales expanded frequency (j^3l) . 

An extended computation of the square of the Weyl tensor (performed as outlined after 
Eq. fHTl) ) leads to C 5C = 0; for the non-oscillatory part of the quadratic term in the 
fluctuations we find instead 



Since the background square Weyl oc t]"^, this computation indicates that the ratio 
SC^/C"^ increases as t]^ in the large scales regime, while it oscillates in the short scales regime, 
with an amplitude that increases as rj^ . 

In the left panel of Figure [6] we present the full numerical evolution of 5C^ /C^ for three 
modes with different momenta. The behavior of the modes that we find here (pure Kasner 
geometry) should be compared with that discussed in the previous Sections, where the initial 
Kasner evolution was followed by an isotropic inflationary stage. In the evolutions shown in 
that case (for instance. Figs. [1] and |4] ) we had isotropization at the time tjso, and we started 
with an initial time of about lO'^tjso. In the present case, the geometry does not undergo 
isotropization. However, the two scale factors are normalized in such a way that they are 
equal to each other at the time t^, cf. eqs. ( I35l) . Therefore, we also choose to = 10~^t* in 
the present case. Also in analogy to the modes shown in Figs. [Hand HI we define k^, to be the 
comoving momentum of the modes which have parametrically the same size as the average 
horizon at the time t^, namely k^^/a^, = 1/ (3t*) (cf. the expressions (136|) ). Moreover, we 
choose ki = kx as in those two Figures. 

Figure [6] confirms that each mode evolves from the non oscillatory large scales regime to 
the oscillatory short scales regime (the transition occurs at later times for modes of smaller 
momenta / larger scales). The time dependence of 5C^ /C"^ shown in the Figure agrees with 
the one obtained analytically, and summarized in Table [H For comparison we also plot in 
the Figure [6] the evolutions of the if+ mode during expansion considered in the Appendix 




, large scales 



, short scales 



(46) 



A(a). 
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- Expanding Kasner - Expanding Kasner 




t / t. t / t. 

FIG. 6: Contribution to 5C'^ /C'^, defined in eq. ([IT]) , from Hy^ (left panel) and if+ (right panel) 
modes with some fixed comoving momentum. The time t* is a reference time at which the two 
scale factors have the same normalization, fc* is the comoving momentum of the modes which have 
parametrically the same size as the two horizons at the time t-t (ki = kx in all the cases shown). 
See the discussions after eqs. (j46]l and (IA4p for more details. 

The results shown in the Figure confirm the instability of the background Kasner solution 
against the GW polarization ■ The growth in the large scales regime (early times) 
agrees with the amplification of the power spectrum shown in Figure [H However, contrary 
to what one would deduce from Figure [H we see that the growth continues also in the 
short scales (late times) regime. We recall that the definition of the power spectrum fl29|) 
contains an arbitrariness in the overall time dependence (since one may have used a different 
combination of the two scale factors as overall normalization). We nonetheless adopted it 
to show the strong scale dependence of the evolution of the power spectrum (which is not 
affected by the overall normalization), and the very different evolution experienced by the 
two GW polarizations (which is also independent of the arbitrary normalization, since Ph^ 
and Ph+ are normalized in the same way). To properly study the instability, one must study 
invariant and unambiguous quantities, such as the (scalar) square of the Weyl tensor which 
is investigated in this Section. 
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2. Hx mode during contraction 

We consider now the contribution to the square of the Weyl tensor from the mode in 
a contracting Kasner geometry. Plugging (1551) into fHUl) . the frequency of the Hy. mode on 
the contracting background has the short and large scales expansions 

2 



(^^) ' short scales 



< - ""^ „ • (47) 



+ i large scales 
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FIG. 7: Contribution to 5C^/C^, defined in eq. ()4ip . from Hx (left panel) and (right panel) 
modes with some fixed comoving momentum. The times and momenta are chosen as in Figure 
[6l While in Figure [6] the background geometry is expanding (from the singularity at t = 0"*"), the 
background geometry is contracting (towards the singularity at t = 0~) in the evolutions shown 
here. 

Once expressed in terms of absolute values of the time, the short and late time asymptotic 
expressions coincide with those of the expanding case, cf. eqs. f H3|) . As in the expanding 
case, the short scales regime occurs asymptotically far from the singularity, while the large 
scales regime occurs asymptotically close to the singularity (notice, however, that a mode 
evolves from the large scales to the short scales regime in the expanding Kasner, while from 
the short scales to the large scales regime in the contracting Kasner background). 

Consequently, also the short and large scales asymptotical solutions are identical, once 
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expressed in terms of t] and rj^, . 



^ 2(-^.) + e"* 2 (-'7,) , short scales 



Ci (^2 kr a/(-?7*) (-??)) + C2 v^^3 (2 /ct a/(-^*) (-^)) ' ^^^g^ scales 

(48) 

However, in contrast to the expanding case, we now keep both solutions for in the large 
scales regime ^. 

For the non-oscillatory part of the Weyl tensor we find 

^^.^^ ^CiQe (-) +h.c.J -^j^^^ , short scales ^ ^^^^ 
M^^X^g'^U)" , large scales 

Notice that the short scales asymptotic evolution agrees with the corresponding one in 
the expanding case. However, in the large scale regime 60"^ is now controlled by the mode 
that was disregarded in the expanding background. 

The asymptotic behaviors fHOl) are confirmed by the fully numerical evolutions shown 
in the left panel of Figure [71 For comparison we also plot the evolutions of the mode 
considered in the Appendix A(b). 



B. Comparison with (in)stability analysis of BKL 

How does the instability of gravitational waves which we demonstrated above fit with the 
classical picture of the universality of the rule of alternation of the Kasner epochs during 
contraction towards a singularity? 

Let us briefly recall this picture. One of the points of the original paper 



was to extend 



the anisotropic homogeneous contracting Kasner solution to a class of generalized Kasner 
solutions, describing more general inhomogeneous anisotropic geometries. It was implicitly 
assumed that there are large-scale, super-horizon inhomogeneities at the scales exceeding 
the (average) Hubble radius. Generic anisotropic solution shall contain eight physically 



different arbitrary functions of the spatial coordinates. However, it was identified by 28 1 



that the homogeneous Kasner solution is unstable against a particular inhomogeneous mode 



^ We recall that this mode was disregarded in the solution (|44|) . since it is a decreasing mode in that case. 
In the contracting case instead this mode dominates at late times. 
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for which / x rot I = (here / is the Kasner axes corresponding to the direction that is 
expanding). In other words, this mode is unstable and is growing with time. Therefore 
the monotonic (stable) inhomogeneous solution can have only seven arbitrary functions of 
the spatial coordinates. Homogeneous Kasner contraction (in Bianchi models with spatial 
curvature) occurs in the stochastic regime of alternation of Kasner exponents 29|. This 
influenced the philosophy of inhomogeneous generalization of the Kasner contraction. Fol- 
lowing [3^, now one can allow all eight arbitrary functions to describe generalized Kasner 
(Bianchi I) geometry. According to jsol , the backreaction of the growing mode of the spatial 
instability alters the Kasner exponents of the patch of the contracting universe in the same 
manner as they were altered in the homogeneous Kasner-oscillating universe. Again, in this 
picture the inhomogeneity scale is larger than the "Hubble" patch, so that locally, along the 
time geodesies, the contraction is asymptotically homogeneous 33| . The instability induced 
by the spatial curvature associated with the inhomogeneous growing mode, and the result- 
ing rotation of the Kasner axes were rigorously studied with the mathematical tools of the 
theory of dynamical systems in 3J, |35|] . 

In the previous Section we found instabilities of both expanding and contracting Kasner 
geometries against gravitational waves. Here for comparison with the BKL analysis we focus 
on the GW instability in the contracting universe. In this background, both polarizations 
Hx and are unstable in the large wavelength limit. All physical wavelengths of the GW 
corresponding to the momenta ki,k2,k3 in different directions (associated with different 
Kasner exponents) sooner or later will leave the Hubble radius ~ ~ t, independently 
of whether they are red- or blue-shifted. The physical frequency of the two modes has the 
structure u;^ ~ + jp- + + f {Ha, Hb), where the function / is of order of H^^- For 
small t this function dominates over the momenta terms; this turns the time dependence 
of the GW amplitude from the oscillating to the non-oscillating regime. The amplitude of 
the non-oscillating GW increases with time in this long-wavelength regime, signaling the 
instability which we described in the previous Section. It turns out (See Appendix [B]) that 
these unstable GW modes in the long-wavelength limit exactly coincide with the unstable 
solution / X rot I = found in [28]. Phrased in another way, the unstable solution of 28] 
which destroys the monotonicity of the homogeneous Kasner contraction is nothing but 
the GW polarizations that are evolving with time into the long-wavelength regime. This 
provides us with a new insight into the origin of the generalized inhomogeneous contracting 
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Kasner solution: small short-wavelength GW will eventually be stretched to become the 
long-wavelength inhomogeneous modes, which are unstable. 

Equipped with the BKL conjecture, one may think that for the contracting universe the 
growth of the unstable GW modes results in an alternation of the Kasner exponents. While 



this conjecture is supported for a contracting universe [2l|, |30(], it is not clear how the GW 
instability evolves for the expanding Kasner background. Indeed, for contracting Kasner 
universe GW leave the Hubble radius and become long-wavelength inhomogeneities, and 
their backreaction can be described by the impact of the spatial curvature at the "local" 
time evolution of a{t) , b{t) , c{t) . For an expanding universe, initially long- wavelength GW 
enter the Hubble radius. Their backreaction can be treated with the pseudo-tensor T^^, = 
Mp < hP^hp^.^ > for the high-frequency GW. 

Another interesting issue is how our unstable modes correspond to the perturbatively 
small variations of the Kasner exponents (giving contribution to the diagonal metric fluc- 
tuations) as well as small variations of the Kasner axes (contribution to off-diagonal metric 
fluctuations). 

Finally, one of the most interesting application of the effect is related to zero vacuum 
fluctuations of gravitons in contracting Kasner geometry, which are unstable. They grow 
and become large scale classical GW inhomogeneities described by the random gaussian 
fleld. One may think that the Kasner axes will be altered and rotated differently in different 
spatial domains of that random fleld. 

We will expand these considerations in a forthcoming investigation. 



VII. AT/T DUE TO THE RESIDUAL CLASSICAL GW FROM PRE-INFLATION 

The growth of the decoupled tensor polarization at large scales can leave an imprint in 
the amplitude and in the polarization of the CMB anisotropics. A characteristic signature is 
a non-diagonal correlation between different multipoles in the expansion of the anisotropies, 
due to initial background anisotropy. Such extended phenomenological study is beyond 
the goals of the present work. However, we want to obtain a crude estimate on the limits 
that such a study would impose on the model, namely on the physical wavelength beyond 
which the statistics of the modes is anisotropic, and on the initial amplitude of the GW 
signal. For this reason, we compute the contribution of the GW mode to the quadrupole 
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of the temperature anisotropies, and we impose that it does not exceed the WMAP value 



C2 ~ 3 • 10-^1 



As we show in appendix [Ul the Ci coefficients of the temperature anisotropies are related 
to the primordial spectrum of the GW by 

^ f "») /: f f t [-"^ « - (^)i ' '») 

where Ig is a "window function", which for a matter dominated universe ^ is given in 
Eq. ( lC12p . The analytic expression for l2{kr]o) can be found in [25]. In this case, the 
function peaks at /c?7o — 3, while it goes to zero as ~ {kfjof' at small argument, and as 
~ cos [krio) I (kriQ)'^ at large argument. The quantity t]q is the present conformal time. For 
a matter dominated universe 

H = — Vo = ^ (51) 

ao 7] ao Ho 

where here ao denotes the present value of the scale factor. The comoving momentum k of 
a mode is related to its present wavelength A by A = 2 vr/ (/c/ao)- Therefore the quantity 



-1 



kr]o = 2n(^^^j (52) 

gives the present ratio between the size of the (particle) horizon and of the mode with 
comoving momentum k . As shown in Figure [21 the power spectrum of the decoupled tensor 
polarization experiences a growth only for k < kiso ■ If the anisotropic stage is followed by a 
prolonged inflationary stage, these large scale modes are inflated to scales much larger than 
the present horizon size, kis.oVo ^ Hq^/Xiso ^ 1 , and the window function suppresses the 
contributions of these modes to the Ce coefficients. 

We now relate kiso Vo to the duration of inflation through the number of e-folds (k) 
between the moment in which the mode with comoving momentum k leaves the horizon 



during inflation, and the end of inflation. This quantity is 22 1 



''' With this assumption we disregard the modes that reentered the horizon during the radiation stage, 
which do not contribute significantly to lowest Ci's, and the period of late accelerated expansion; this is 
consistent with the present approximate computation. 
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For a matter dominated universe aoHo = 2/r]o, cf. eq. ( 15T|) : Vq is the energy density in 
the universe at the moment in which that mode left the horizon (for the potential we are 
discussing, cf. Eq. (|2]), this quantity is nearly independent of k); finally, the quantity A is 



sensitive to the details of reheating (see 23|] for a discussion). Imposing that all the modes 
within the horizon today were sub-horizon at the beginning of inflation gives the minimal 
number of e-folds given in Eq. ([1]). We define 

AN^N (fci,o) - iV^in = -In = -In (^^) (54) 

Since kiso is the comoving momentum of the modes leaving the horizon when the universe 
becomes isotropic, AA^ is the difference in e-folds between the duration of the isotropic 
inflationary stage and the minimal duration ([1]). As we discussed, a large AA^ results in a 
suppressed effect of the amplified GW modes on the observed CMB anisotropies. 

We compute only the contribution of to the C2 coefficient in Eq. (l50ll . Figures [2] and 
[3] give the power spectrum of this mode at the end of inflation normalized by the initial 
power spectrum Pq. This quantity is related to the initial condition at the start of the 
anisotropic stage, and it was left unspecified. For definiteness, we parametrize it with an 
overall amplitude and with a power law dependence on the scale 

Po = Aox{7]okr , (55) 

where ^0 is the initial power at a scale which parametrically corresponds to the present 
horizon size. We numerically perform the angular integral in (1501) . to find 



\hsoJ J-l [ (^)' , k > fci,o 

where the proportionality to \ho\ is related to the growth of the power spectrum during the 
anisotropic era, as explained at the end of Section HVl The dependence on k is due to the 
fact that the angular integral is dominated by modes of small angle {ki ^ k^), cf. Figure 
131 where the growth of the spectrum is flat for k kiso, and scales as ~ k~^ at larger k . 

The function P is shown in Figure [H normalized by |/io| and by the initial power spectrum 
Pq . The dependences on \ho\ and on k given in eq. ( !56|) are manifest in the Figure. 

We can now perform the final integral over momenta in (l50l) . Taking into account the 
large and small argument dependence of I2 and that of P given in eq. (l56l) . we see that the 
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FIG. 8: Angular integral of the power spectrum, defined in Eq. (j56p . We show the result normalized 
by the initial power spectrum, starting from different values of |/io| (given in units of -v/Vo/^oi 
larger values of this quantity corresponds to a longer anisotropic stage), and rescaling it by \ho\, 
to explicitly show the linear dependence of the power spectrum on this quantity. 

integrand behaves as ~ k"''^ at small k, and as ~ at large k. Therefore, the integral 
converges for a wide range for the initial slope defined in Eq. ( l55l) . namely —4 < n < 7 . 

After performing the integral, we impose that the resulting C2 does not exceed the WMAP 
value C2 — 3 ■ 10^^^ [24]. This results in an upper limit on ^ol^ol (the initial overall 
amplitude, times the growth during the anisotropic phase) for any value of fciso and n . We 
show this in Figure [9l As explained after Eq. (152|) , the limit weakens at small values of 
hsoVo, corresponding to a longer duration of the inflationary stage. 

Finally, we show in Figure [ID] the first few C/s, normalized by the initial amplitude of 
the GW and the length of the anisotropic era (controlled by |/io|)) for a few values of the 
slope fl55|) . As expected, for moderate slopes, the spectrum of Ci decreases with i. This 
is due to the fact that larger angular scales (small i) are affected by the modes with lower 
momenta, which grow more during the anisotropic stage. The three spectra shown in the 
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FIG. 9: Limit on the initial amplitude of the power spectrum, times its growth during the 
anisotropic stage. The quantity AA^, defined in eq. ()54|) . is the number of e-folds of isotropic 
inflation minus the minimum usually required to homogenize the universe, cf. eq. ([T|). For a longer 
duration of the inflationary stage, the modes experiencing the growth are inflated to larger scales 
than the present horizon, and the resulting effect is suppressed (weaker limit), and n, are, 
respectively, the amplitude and the slope of the initial power spectrum, see eq. (f55|) . while \ho\ is 
proportional to the duration of the anisotropic stage. 

figure can be fitted by a single power-law (with an accuracy up to about 10%): 



80^3-8 , n = -2 
248 r , n = 
1430 r°-^ , n = 2 



(57) 



We see that the angular spectrum roughly decreases with £ as ~ l/i^, where the 
exponent p depends on the power-spectrum of the classical GW at the initial hypersurface. 
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FIG. 10: First Cg coefficients generated by the GW, normalized by tfie initial amplitude and the 
duration of the anisotropic era (controlled by |/io|)- The tree curves corresponds to three different 
slopes of the initial power in the GW, eq. ([55|) . 



VIII. DISCUSSION 

We found a new effect of instability of the gravitational waves in an expanding and con- 
tracting Kasner geometries. We demonstrated the effect for a particular choice of the Kasner 
exponents (—1/3, 2/3, 2/3), but we expect it is rather generic. This particular choice allows 
to simplify the description of GW in such a way that their wave equations is manifestly 
depending on the momenta. For the contracting Kasner geometry, we found that our unsta- 
ble GW mode is identical to the unstable large-scale inhomogeneous mode first identified in 



28|-|30| (for arbitrary Kasner exponents (pi,P2,P3))- Backreaction of this unstable mode is 



conjectured to alter the Kasner exponents, differently in different spatial patches, depending 
on the spatial profile of the initial GW. 

The Kasner geometry is a rather universal asymptotic solution for many interesting sit- 
uations: it describes approach towards black holes or cosmological singularities (including 
higher dimensional and supersymmetric cases [36]), and it describes generic anisotropic ex- 
pansion from singularity prior to inflation. 
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All of these situations are inevitably accompanied by quantum fluctuations of the gravi- 
tons), and, possibly, also by classical gravitational waves. There is a long list of questions 
arising in connection with a new effect of the gravity waves instability in anisotropic geom- 
etry. We have to investigate how the instability growth depends on the Kasner exponents 
PijP2,P3- It will be interesting to understand if the instability of the classical GW also 
results in the instability of the graviton zero vacuum fluctuations. It will also be worth to 
understand the impact of the GW instability on the structure of the singularity inside the 
black hole and cosmological singularity. It is also interesting to investigate the backreaction 
of the GW instability on the contracting and expanding anisotropic geometries depending 
on the initial GW proflle. 

In this paper we considered the effect of GW instability in the context of the anisotropic 
pre-inflationary stage. The transition from Kasner expansion to inflation terminates the 
effect of GW instability but leaves classical GW signal as the initial conditions for inflation. 
If inflation does not last very long, then the residual GW can contribute to the CMB 
temperature anisotropies. Since GW polarizations and the scalar mode of cosmological 
fluctuations behave differently during anisotropic pre-inflation, we can consider only the 
leading contribution, namely, mode of the GW polarization. In this paper we calculated 
its contribution to the total AT/T anisotropy angular power spectrum. The angular power 
spectrum of the signal decreases with £ in a power-law manner and depends on the initial 
spectrum of the classical GW fluctuations, Ce ~ l/i^, where the exponent p depends on the 
power-spectrum of the classical GW at the initial hypersurface. It is interesting to note that 
this result is qualitatively similar to the result of 37|, where the impact of GW from the 
our-universe-bubble nucleation was considered. 

The signal from is rather anisotropic, and there is an interesting question about the 
anisotropy of its multipole structure {aimdl,^,), which is intriguing in connection with an 
apparent alignment of the low multipoles of AT/T. While such an analysis is beyond the 
aims of the present work, we have estimated the initial conditions (initial amplitude of the 
GW, versus the duration of inflation) which can lead to potentially observable effects. We 
leave a more extended analysis to future investigation. 



Acknowledgements 



34 



We thank J.R. Bond, CR. Contaldi, T. Damour, I. Khalatnikov, A. Linde, C. Pitrou, 
M. Sasaki, A. Starobinsky, J. P. Uzan and J. Weinwright for useful discussions and cor- 
respondence. The work of A. E.G. and M.P. was partially supported by the DOE grant 
DE-FG02-94ER-40823. LK was supported by NSERC and CIFAR. 



APPENDIX A: CONTRIBUTION OF THE H+ MODE TO THE WEYL 
INVARIANT 

In this Appendix we compute the contribution of the if+ mode to the square of the 
Weyl tensor for a pure Kasner background. We treat the expanding and contracting cases 
separately. This presentation follows the analogous one for the ifx mode placed in Subsection 
IVI Al of the main text. 



a. mode during expansion 



Substituting fl35|) into fl40|) we have the large and short scale expansions 

large scales 
\ — \ ■. short scales 

Consequently, we have the asymptotic solutions 




(Al) 



Ci^Jo (2 kr y/rpq) , large scales 

r c -.^ , , • ^^^^ 

■^te 2»7. + i^e , short scales 

where again in the early time solution we have disregarded a decaying mode that would 
diverge at the initial singularity. Notice that the late time expression is the same as for the 
Hx mode (since the late time expansions of uj\ and uj\ coincide). 
For the linear part of the Weyl tensor we find 

, ^^^^li^^ , large scales 



Cie +C2e 2,, ) , short scales 

The non-oscillatory part of the 2nd order term in the Weyl square is 

Af 2 al rf* ' 1^^§6 scales 




Ci C2 e "* + h.c. ) j2^{iA^5 , short scales 



(A4) 
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This behavior is confirmed by the fully numerical solutions shown in the right panel of 
Figure [61 It is interesting to compare the contribution of the mode and that of the iJ+ 
mode. In the large scale regime only the contribution of the i^x niode grows with respect 
to the background. This is consistent with the amplification of this polarization in the 
Kasner era that we have found in Section IIV[ In the small scales regime, the contribution 
of both modes is instead growing with respect to the background. Although the two modes 
evolve identically in this regime (since they satisfy the same asymptotic equation), the time 
evolution of the corresponding 5(7^ term is different, since the two modes enter differently 
in the Weyl tensor. 



Hj^ mode during contraction 



Plugging fl35l) into fHOl) we have the early / late time expansions 

2 



kr ( I , short scales 

ujI^I ^ v-"*; . (A5) 

i T(~n)' ^^T^ , large scales 
The similarities and differences between the if+ mode in the expanding and contracting 
cases are identical to those discussed in the previous Subsection for the mode. We have 
the asymptotic solutions 



H. 



2{- 



v^) -\- e 2^-^*) , short scales 

-1,/— rj * 



Ci -^0 (2 kT V(-?7*) (-??)) + C2 ^/^Yo (2 kr ^/ {-rj*) , large scales 

(A6) 

This leads to the following asymptotic evolution for the linear term in the square of the 
Weyl tensor 



Ci e ^c-..) + C2 e" ^'Y^^V, \l , short scales 

Alp a^i (—11) 



SttCi + 6^2 (in (kT ^A^V*)i^) + 7 - 



>+h.c. 



X ^^yy^ .6 , large scales 



(A7) 

and for the non-oscillatory part in the quadratic term of the square of the Weyl tensor 

CiGne "J* + h.c. f,,A, .5 , short scales 
SC''^{ \ _\ ^ J ^^^^(-'^.r _ (A8) 

M2^!2a6 4 ' ^^^^^ ^^^^^^ 
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This behavior is confirmed by the fully numerical evolutions shown in the right panel of 
Figure [3 



APPENDIX B: INSTABILITY OF KASNER SOLUTION: KL VS GW 



Our analysis of the squared Weyl invariant shows that the Kasner background is unstable 
due to the amplification of the GW perturbations and ifx • This is manifest by the growth 
of 5C/C and bC'^ jC"^ summarized in table [Tl For the contracting background we found that 
the contribution to bC"^ jC"^ of the polarizations Hj^ and iJx grow, respectively, as 
and as asymptotically close to the singularity. Conversely, in the expanding case, the 
Weyl invariant is regular near the singularity. The reason for this discrepancy is that, in 
the expanding case, we disregarded the Neumann term Yq and Y3 (for ifx) iu the 

solutions of the evolution equations (1391) for the two modes. However these two modes are 
generally excited as the contracting Kasner background solution approaches the singularity. 

The purpose of this Appendix is to show explicitly that this instability is not inconsistent 
with the lore of alternation of Kasner exponents l29| . On the contrary, the instability in 
the contracting case was pointed out already in [28j, where it was shown that the Kasner 
background is stable provided one physical condition is imposed. As we now show, the 
physical condition imposed by |28|] precisely eliminates the unstable Yq and Y3 solutions. 



We start by quoting the results of [28| where the authors considered a perturbed vacuum 



Kasner solution in synchronous gauge and did a stability analysis near the singularity. In 



28|, 



the main text, we describe the contracting Kasner background using negative time. In 

instead, positive time was used, and the approach to the singularity was studied in the 

i 

t ^ 0^ limit. In this Appendix, we adopt the time convention of [2^]. Moreover, we specify 
the analysis of 28|] to the axisymmetric = —1/3, P2 = Ps = 2/3 case. For this choice, the 



solution obtained in 



28j for t ^ 0+ reads 



hi 

h23 



Ai + Bi \nt 



hi = A2 + B2 Int 



^-2/3 



9 



Cl2--fc3 {k2Ci3-hCu)t'/' 

Ci3--^k2 {hC,2-k2Cr3)t"'' 



t'/'' {C23 + C; 



32 J 



hi = A3 + 53 Int 



21 



31 ; 



(Bl) 
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with constants Ai,Bi, Cij satisfying the constraints 



B^ + B2 + Bs = -^B, + ^ (52 + 53) = , 
h (B, -A2-A,) + 2 {k2 C21 + h C31) = , 
k2{B2 + A^)-2hCu = 0, 

ksiB, + A,)-2hCis = 0. (B2) 



As a result of the synchronous gauge, the solutions fIBip still have freedom due to gauge 
artifacts; the corresponding gauge contributions to the perturbations are given in eqs (F3) 



of 



281]. Using this freedom, the Authors of 



28j set C12 = C31 = C23 + C32 = 0. « After this 



gauge fixing, one sees that the mode hi^ diverges at the singularity. For this reason, ref. 



28l ] concluded that the Kasner background is stable only provided the condition C13 = is 
imposed. We stress that setting C13 is not a gauge choice, but rather the suppression of a 
physical degree of freedom. This physical constraint, along with the relations (]B2p . implies 
that B2 = 0. 

We now turn back to our analysis. For t — > 0+, eqs. (l39l) admit the solutions 

= t^'^ [E,U?,kTt^l^) + E2Y^{?,kTt^l^)\ . (B3) 

The modes Jq and J3 are regular at the singularity, whereas Yq and Y^ diverge, thus resulting 
in the instability. We now show that the conditions C12 = C13 = B2 in I28] correspond to 
having D2 = E2 = Q m. our case. To do such a comparison, we need to change from the 



gauge chosen in the main text to the synchronous gauge adopted in [28[. We do so through 
a general infinitesimal transformation — > x'^ + , with E,^ = {^^, di C,^, d^^ + where 
is a transverse 2d vector. From the transformed metric 6gfj,u 6g^i, — C^g^J , we find that 



They are in fact setting C23 = in the general case. However, the time dependences of both terms in h23 
are the equal to each other for the axisymmetric p2 = Pa case; therefore, in this case, the removal of the 
gauge artifact corresponds to C23 + C32 = 0. 
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the parameters ^'^ which take our gauge to the synchronous coordinates are 

e = — fdt'm +—f, 



e = 
e = 



dt 



, x{t') 

a{t') 



+ 



■ dt' 
* dt' 



h\t') 



dt"^{t" 



dt" <^{t" 



^ J a^t')^^ 



(B4) 



where /, /* are integration constants and fcj/* = 0. The perturbations $, x, S and 

Bi are the modes in the decomposition given in the main text. Using their relations to 
the canonical modes if+ and H^, along with the early time solutions ( 1B3I) . we calculate 
the metric perturbations in synchronous gauge. The relevant components for the present 
discussion are 

12 



'''' - ~27^^wm; 



2/3 I ^ 1.2 , ^ . 4 ,2/3 

8 ^ 64 ^ 



fk^ k2t + ki k2 + 0(t4/3) 



il3 



hi 



16iV2k2 
27 71 k^Mp 



E2 {t-^'' + ^ + ^ 4 t^''") -ffk,k,t + k^ k, f + 0{t'l^) 



V2 



3nM„ 



371 Di + 2 - 2) D2 + 6 D2 In 



3ki 



2V2 



277rk^Mp 



k 

32 i kl k2 -^E2 - 9 Da 4 {^1 - kl) 
kr 



Int 



-2 /o (^3 kl + + 2 ^2 /2 + fc2 /) + Oit'/') . 



(B5) 



In eqs fIBip . setting C12 = (gauge choice) and C13 = (physical choice) eliminates the 
terms up to 0{t^^^) in hi2 and hi^. On the other hand, from our solution flB5l) . we see that 
these terms in hu vanish only if -E2 = 0. Furthermore, the choice C12 = C13 = implies that 
B2 = 0, now eliminating the log term in h^, which corresponds to 1^2 = in our solutions. 
In other words, the stability conditions derived in [28j correspond to removing the Neumann 
functions in the early time solutions ( 1B3I) of both and Hx ■ 



APPENDIX C: EXPLICIT COMPUTATION OF C^ ™™/ 



Here, we prove the relation (150|) of the main text. This expression gives the Ci coefficients 
of the temperature anisotropies in terms of the value of the power spectrum of GW at the end 
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of inflation. Therefore, this computation is performed in an isotropic universe. The effects 
of the anisotropic stage are encoded in the power spectra that, on large scales, depend also 
on the orientation of the modes. 

We start by decomposing the temperature perturbations in direction p measured by an 
observer at Xq and at conformal time rjo, into spherical harmonics 

(p) = X] YimiP) ■ 



T 



(CI) 



Assuming that the perturbations are Gaussian, their statistics are completely specified by 
the second order correlations, which can be written by inverting (IC1|) 



I dnpdnp, ^^(p) ^(p')^ Y^^{j>)Y^.m'{v') . (C2) 



Assuming instantaneous recombination, we can use the well-known solution to the first order 
Boltzmann equation for gravity waves: 

1 

~ ~2 



dv - 

'?LSS ^ 



(C3) 



where hij = dgij/a^. We expand the gravitational perturbations in terms of plane waves as 



d^k 



A=+,x 



(2 7r 



,3/2 



ik-x (A) 

ij 



(k) h^'Hk,^)ai'^ 



(C4) 



where (A) designates the two polarizations (+) and (x), e-^'*(k) is the normalized polarization 
tensor satisfying e\j\\<.)e^^y^(k) = S^^y^, and a^^' are Gaussian random variables with unit 
dispersion, ie. {o-k^ ^ {x')) — ^'''^^(k — h') 6^^yy Using the solution ( ]C3I) along with the 
decomposition (lC4p in the expression (]C2|) . we get 

d^k 

4^./ (2 7r)3 



(C5) 



Cff'r 



dQp dVtpi Y;,^(p)Y£/m'(p 



where we noted that x = xo + p(77 — r/o) describes the geodesies of the photons with direction 
p, observed at (xo,r7o)- 

The contraction e\^'' (k) fp can be calculated in terms of the angles of p defined with 



respect to a coordinate system with a 2;-axis coinciding with k, as 



26| 



e^{k)p^f = 5l^J^ cos(20^.) + S[:\ sin(20, 



c(A) 



(A) 



sin 9ri 



(C6) 
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To take the angular integral over dQp (dflpr), we first rotate the spherical harmonics to the 
angular basis (p'^) with 

Yirnip) = J2 ^-M(k) Ymip^) , (C7) 

M 

where -D^\^ are the Wigner coefficients satisfying '^M^Mm^Mm' — ^mm'- Finally, the 
angular integrals can be evaluated as [261] 



'{ 2i+l){i + 2)\Je+i/2 [k ivo-v)] 
2(^-2)! [k{vo-ri)f' 

(C8) 

Ignoring the effect on the (small scale) modes which enter the horizon during the radiation 
dominated era, the solutions for t he g ravity waves during matter domination can be written 
in terms of the primordial one as 25], 



Using the integral (ICSp along with the above solution yields 

_ V- , /(2£ + l)(2£' + l)(£ + 2)!(f + 2)! 



(C9) 



-I] 



16 y (£-2)!(f-2)! 

X \{^M ,2^M' ,2 + -2^M> ,^2) (<^A +^a) ^ M -2^1' ,2 + ^ M ,2^ M> -2) (^A ~^a)] 

r T /r f ^ ^^'^^'^^^ ^^^^^^ ' 

where ^ = cosOk, the power spectrum is defined through 

P„^„).«M)P, ,,,,, 

which coincides with the definitions given in the main text, cf. eqs. fl29p and (1331) once the 
universe has isotropized. Finally, the time integral is 

hikr]o)= dy—jj^— (C12) 

To be able to compare with observations, we average the diagonal {^ = m = m') part of 
the correlator over all m and obtain 



/r f r ^ ^^''^^''^ ^ ^''^ ^^^^ ' ^^^^^ 
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The final result (1C13P is valid for any model of pre-inflation which results in an anisotropic 
primordial tensor power spectrum. 

For standard (isotropic) inflation, the power spectrum is the same for both polarizations 



and we recover the standard result 



25|. 
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